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1. Introduction 

In this work we compare the spectrum of the inverse decay lengths of Polyakov loop correlators 
in the {3+V}d SU(3) gauge theory in the deconfined phase near the transition with the spectrum of 
massive excitations of the 3<i 3-state Potts model in the broken phase near the transition at zero 
magnetic field, which were determined in Ref. ||l|]. The aim of the work is to verify if and to what 
extent the Svetitsky-Yaffe conjecture also holds for theories which undergo a weakly first order 
phase transition, using mass ratios as a probe. In particular we focus on the low-lying masses 
in two different sectors of parity and orbital angular momentum, 0+ and 2+. We expect that, if 
universality would apply in strict sense, these spectra should exhibit the same pattern, as suggested 
by several numerical determinations in the 2x1 Ising class [^, 0, ^ |6|]. 

We extend our numerical analysis to temperatures far away from the transition temperature T, 
in order to look for possible "scaling" of the fundamental masses with temperature. Moreover, we 
consider also the screening masses resulting from correlators of the real parts and of the imaginary 
parts of the Polyakov loop. These determinations can represent useful benchmarks for effective 
models of the high-temperature phase of SU(3), such as those based on mean-field theories of the 
Polyakov loop, suggested by R. Pisarski [^. 



2. Screening masses from Polyakov loop correlators 



Screening masses are defined as the inverse decay lengths of the Yukawa-like potential be- 
tween two static sources. They are generally determined through the correlation of suitable opera- 
tors. In our case correlations are between operators with different spatial separation. The general 
large distance behavior for the correlation function G(|zi — Z2I). in an infinite lattice, is: 



G(ki-Z2l)=I«„e-'""l^'-^^l 



(2.1) 



where tuq is the fundamental mass, while m\, m2, ... are higher masses with the same angular 
momentum and parity quantum numbers of the fundamental mass. On a periodic lattice the above 
equation must be modified by the inclusion of the so called "echo" term: 



G{\Zl -Z2I) = 52^n 



^-m„|zi-z2| _|_g-mn(A';-|zi-Z2|) 



(2.2) 



The fundamental mass in a definite channel can be extracted from wall-wall correlators by looking 
for a plateau of the effective mass at large distances, 

G{z-l) 



meff (z) = In ■ 



G(z) 



In the 0+ -channel, the connected wall-wall correlator in the z-direction is defined as 



G(Izi-Z2l)=Re(P(zi)P(z2)^)-|(P 



where 



Hz) 



L L P{nxa,nya,z) 



(2.3) 



(2.4) 



(2.5) 
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represents the Polyakov loop averaged over the x3'-plane at a given z- ^ The wall average implies 
the projection at zero momentum in the xj^-plane. 



For the 2+-channel, we used the variational method [Q, (for more details, see [|T^] and 
references therein.) 

Our choice of wall-averaged operators in the 2+ -channel is inspired by Ref. [|ri|] and reads 



P{nxa + na,nya,z) — P{nxa,nya + na,z) 



(2.6) 



In most cases we have taken 8 operators, corresponding to different values of n, with the largest n 
almost reaching the spatial lattice size Njc 

We consider also correlators of the (wall-averaged) real and imaginary parts of the Polyakov 
loop, defined as 

GrUzi -Z2\) = (ReP(zi)ReP(z2)) - (ReP(zi))(ReP(z2)) , (2.7) 
Giilzi -Z2\) = (ImP(zi)ImP(z2)) ■ (2.8) 

The corresponding screening masses, mr and mj, can be extracted in the same way as for the 
0+ mass. We have studied the ratio mj/rriR over a wide interval of temperatures above the tran- 
sition temperature Tf of (3-1-1 )J SU(3) and seen how it compares with the prediction from high- 
temperature perturbation theory, according to which it should be equal to 3/2 [12, |l^, and with the 



prediction from the mean-field Polyakov loop model of Ref. according to which it should be 
equal to 3 in the transition region. The interplay between the two regimes should delimit the region 
where mean-field Polyakov loop models should be effective. 

3. Numerical results 

We used the Wilson lattice action and generated Monte Carlo configurations by a combination 



of the modified Metropolis algorithm [ ]15P with over-relaxation on SU(2) subgroups []160. The 
error analysis was performed by the jackknife method over bins at different blocking levels. We 
performed our simulations on a 16^ x4 lattice, for which j3, = 5.6908(2) over an interval of B 
values ranging from 5.69 to 9.0. 

Screening masses are determined from the plateau of mgff (z) as a function of the wall separa- 
tion z- In each case, the plateau mass is taken as the effective mass (with its error) belonging to the 
plateau and having the minimal uncertainty. We define plateau the largest set of consecutive data 
points, consistent with each other within la. This procedure is more conservative than identifying 
the plateau mass and its error as the results of a fit with a constant on the effective masses /Mgff(z), 
for large enough z- 

Just above the critical value jS^ we find a large correlation length, which is not of physical 



relevance. It is instead a genuine finite size effect [18] related to tunneling between degenerate 
vacua. This effect disappears by going to larger lattice volumes or moving away from j8f in the 
deconfined phase. Tunneling can occur between the symmetric and the broken phase, and between 



Here and in the following, Nj (i = x,y,z) is the number of lattice sites in the /-direction. 
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Figure 1: Effective mass in the 0+ (left) and the 2+ (right) channel as a function of the separation between 
walls on the {x^y) plane at j3 =5.75. 



the three degenerate vacua of the deconfined phase. When tunneling is active, the correlation 



function has the following expression 1 18]: 

G{\zi 



■Z2\ 



-mo\zi-z2\ 



+ boe 



-m,\z\-Z2\ 



(3.1) 



where is the inverse of the tunneling correlation length and is generally much smaller than the 
fundamental mass /mq and therefore behaves as a constant additive term in the correlation function. ^ 
The dependence on lUt in the correlation function can be removed by extracting the effective mass 
by use of the combination 

G{z)-G{z+\) 



'"eff (z) = In ■ 



(3.2) 



G{z+l)-G{z + l) ' 

A typical example of the behavior of the effective mass with z is shown in Fig. [T] for the 0+ and the 
2+ channels. In Figs. || and |3| we show the behavior with p of mo+, m2+, and ni/. 

We observe from that mo+ and niR are consistent within statistical errors, this indicating that 
the Polyakov loop correlation is dominated by the correlation between the real parts. We can see 
that the fundamental mass in the 0+ channel, as well as Ar, becomes much smaller than 1 at jSf, 
as expected for a weakly first order phase transition. In the cases of /mo+ and of niR we have made 
some determinations below j8f (see Figs. ^ and |3|). It turns out that masses in lattice units take their 
minimum value just at j8,, where there is a "cusp" in the j8 -dependence. Such a behaviour was 



observed also by the authors of Ref. []19|], whose results, when the comparison is possible, agree 
with ours. We have also looked for a scaling law for the fundamental mass in the 0+ channel, but 
with the understanding that any second-order-like scaling law, when applied to the region near a 
first order phase transition, should be taken as an effective description, which cannot hold too close 
to the transition point. With this spirit, we have compared our data with the scaling law 



(3.3) 



^In ( [3.l[ ) we have taken into account only the lowest masses in the spectrum and, for brevity, omitted to write the 
"echo" terms. 
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Figure 2: Screening mass in the 0+ channel (left) and in the 2+ channel (right) vs. j3 . 
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Figure 3: Screening masses niR (left) and m/ (right) vs. j3. 



where mo+(j3i) and mo+{^2) are the fundamental masses in the 0+ channel at j8i and ^2, respec- 
tively. We have considered several choices of jSi and found that for each of them there is a wide 



"window" of j8 values above j8f where the scaling law (3.3) works, with a "dynamical" exponent 
V (see Ref. [jT^ for a details). For j8i=5.72 we have calculated also the ;^^/d.o.f. when v is put 
exactly equal to 1/3 (suggested in Ref. pO| ] to apply to the standard correlation function), getting 
X^ld.o.f.=0.15 in the window from j3 = 5.715 to j3 = 5.78. In Fig. | we show, for this choice of jSi , 
the comparison between data and the "scaling" function with v set equal to 1/3. 



Then, we have considered the j3 -dependence of the ratio m2+/nio+, shown in Fig. |5[ We have 
found that this ratio can be interpolated with a constant in the interval from j8, to j3 = 5.77. This 
constant turned out to be 3. 172(65), with a ;^^/d.o.f equal to 1.085. In the fit we excluded the point 
at j8=5.695, for which the determination of m2+ is probably to be rejected. If the point at j3=5.695 
is included, the constant becomes 3.214(64) with ;^^/d.o.f =2.21. The fact that the ratio m2+/mo+ 
is compatible with a constant in the mentioned interval suggests that m.2+ scales similarly to mo+ 
near the transition. This constant turns out to be larger than the ratio between the lowest massive 
excitations in the same channels in the broken phase of the 3(i 3-state Potts model, which was 
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Figure 5: (Left) Ratio OT2+ /mQ+ as a function of j3 in the deconfined phase. The three upper horizontal lines 
represent the constant (with its error) which fits the data (see the text for details); the three lower horizontal 
lines represent the corresponding mass ratio (with its error) found in the 3(i 3-state Potts model [|lj|. (Right) 
Ratio mj/ nig as a function of j3 in the deconfined phase. The vertical line corresponds to the critical j3 value. 



determined in Ref. ^ to be 2.43(10). 

We have calculated the ratio nij/mn for j3 ranging from 5.695 up to 9.0. We observe from 
the right panel of Fig. || that this ratio is compatible with 3/2 at the largest j8 values considered, 
in agreement with the high-temperature perturbation theory. Then, when the temperature is low- 
ered towards the transition, this ratio goes up to a value compatible with 3, in agreement with the 
Polyakov loop model of Ref. [p^, which contains only quadratic, cubic and quartic powers of the 
Polyakov loop, i.e. the minimum number of terms required in order to be compatible with a first 
order phase transition. The same trend has been observed also in Ref. [19]. 



4. Conclusions and outlook 



In this work we have studied in the (3-i-l)<i SU(3) pure gauge theory above the deconfinement 
transition the lowest masses in the 0+ and the 2+ channels of angular momentum and parity and the 
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screening masses resulting from the correlation between the real parts and between the imaginary 
parts of the Polyakov loop. The behavior of the ratio between the masses in the 0^ and the 2"*" 
channels with the temperature suggests that they have a common scaling above the transition tem- 
perature. This ratio turns to be ~30% larger than the ratio of the lowest massive excitations in the 
same channels of the 3J 3-state Potts model in the broken phase. This can be taken as an estimate 
of the level of approximation by which the Svetitsky-Yaffe conjecture, vaUd in strict sense only for 
continuous phase transitions, can play some role also for (3+l)J SU(3) at finite temperature. 

The dependence on the temperature of the ratio between the screening masses from the corre- 
lation between the real parts and between the imaginary parts of the Polyakov loop shows a nice 
interplay between the high-temperature regime, where perturbation theory should work, and the 
transition regime, where mean-field effective Polyakov loop models could apply. 
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